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We show that the single, non-photonic electron nuclear modification factor R\a is affected by 
the thermal enhancement of the heavy-baryon to heavy-meson ratio in relativistic heavy-ion colli- 
sions with respect to proton-proton collisions. We make use of the dynamical quark recombination 
model to compute such ratio and show that this produces a sizable suppression factor for R'aa 
intermediate transverse momenta. We argue that such suppression factor needs to be considered, 
in addition to the energy loss contribution, in calculations of R'aa- 
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I. INTRODUCTION 



The suppression of single, non-photonic electrons at 
RHIC [ll, Q is usually attributed to heavy-quark energy 
losses. However, calculations that successfully describe 
the nuclear modification factor of hadrons fail to de- 
scribe the single, non-photonic electron nuclear modifi- 
cation factor R\a S 13; 111 ■ This has prompted a great 
deal of effort aimed to better describe the heavy-quark 
energy loss mechanisms to include not only the radiative 
part [y, 0, [1| but also the coUisional Q and the medium 
dynamical properties to compute the radiative piece . 
As a result, although some improvement in the descrip- 
tion of the nuclear modification factor has been gained, 
it is not yet clear whether the anomalous suppression can 
be completely attributed to energy losses. 

Working along a complementary approach to describe 
the n on-p hotonic electron yield at RHIC, it has been ar- 
gued [ill , [l^ that under the assumption of an enhance- 
ment in the heavy-quark baryon to meson ratio, analo- 
gous to the case of the proton to pion and the A to kaon 
ratios in Au-|-Au collisions [Tsl. [Tsl. [l6t . it is possible to 
achieve a larger suppression of the nuclear modification 
factor. The rationale behind the idea is that heavy-quark 
mesons have a larger branching ratio to decay inclusively 
into electrons as compared to heavy-quark baryons, and 
therefore, when the former are less copiously produced 
in a heavy- ion environment, the nuclear modification fac- 
tor decreases, even in the absence of heavy quark energy 
losses in the plasma. 

In order to give a qualitative argument that shows how 
an enhancement in the heavy-quark baryon to meson ra- 
tio can suppress the single, non-photonic electron nuclear 
modification factor, let us look at the px integrated R\a 
and to consider that the heavy hadrons are only those 



containing a single charm, 
1 
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where (rip) is the average number of participants in the 
collision for a given ccntrality class, refers to 

the number of x-particles produced in A -f A (p + p) 
collisions and B^^^ is the branching ratio for the in- 
clusive decay of x-particles into electrons. Let us bring 
Eq. dll) into a form that contains the corresponding pT 
integrated nuclear modification factor for particles con- 
taining charm. We write 
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Let us introduce the shorthand notation 
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where C represents the enhancement factor for the ratio 
of charm baryons to mesons in A -|- A as compared to p -|- 
p collisions and x is the charm baryon to meson relative 
branching ratios for their corresponding inclusive decays 
into electrons. With these definitions, and after rewriting 
the factor N^^^/N^ in the form 
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FIG. 1: (Color online) pr integrated T2_a as a function of x, 
the ratio of branching ratios for charmed baryons and mesons 
to decay inclusively into electrons. Notice that for x < 1, 
TXa < 1 when Ca -the ratio of charm baryons to mesons in 
A + A- is larger than a -the ratio of charm baryons to mesons 
in p + p-. 



we can express Eq. ^ as 
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When not integrated over transverse momentum, 
the factor 1/ {n^) [{N^^ + N^^)/{Ng + N^^)] reprc- 
sents the nuclear modification factor for particles with 
charm. Let us not assume any particular value for this 
factor and instead concentrate in the other one in Eq. ([5]), 
which can be written as 
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where a = Np^/N^p. The above quantity is plotted in 



Notice that the function T^^^'°' is less than 1 



Fig. [T] as a function of x for different combinations of Ca 
and a. 

when X < 1 provided that Ca > a 

In this work, we want to quantitatively address the 
question of whether the enhancement factor C times a - 
namely, the hcavy-baryon to heavy-meson ratio in Au + 



Au collisions- can indeed be larger than a -namely, the 
heavy-baryon to heavy-meson ratio in p -I- p coUisions- 
and if so, how this affects the behavior of the factor 
as a function of p^. For these purposes, we use a dynami- 
cal recombination scenario that accounts for the fact that 
the probability to form baryons and mesons can depend 
on a different way on the evolving density during the 
collision. 

The work is organized as follows: After presenting a 
brief introduction to the dynamical quark recombination 
model in Sec. [Ill we proceed in Sec. IIIII to compute the 
probabilities to form mesons and baryons containing a 
heavy quark in a relativistic heavy-ion collision environ- 
ment. In Sec. IIVI wc use these probabilities to write ex- 
pressions for the meson and baryon transverse momen- 
tum distributions. In Sec. |V]wc compute such distribu- 
tions as well as the baryon to meson ratio. We convolute 
such ratio with the branching ratios of charmed baryons 
and mesons to decay into electrons to obtain the pt un- 
integrated function and show that this can be in- 
deed less than 1. Finally we summarize and conclude in 
SedVH 



II. DYNAMICAL QUARK RECOMBINATION 

Recall that hadronization is not an instantaneous pro- 
cess. In fact, lattice calculations [Ol show that the phase 
transition from a deconfined state of quarks and gluons 
to a hadron gas is, as a function of temperature, not 
sharp. Working along this line of thought, it has recently 
been shown [l8| that the features of the proton to pion 
ratio can be well described by means of the so called 
dynamical quark recombination model (DQRM) that in- 
corporates how the probability to recombine quarks into 
mesons and baryons depends on density and tempera- 
ture. Other approaches toward a dynamical description 
of recombination have been recently formulated [l9| . 

The upshot of the DQRM is that the density evolving 
probability differs for hadrons made up by two and three 
constituents with the same mass, that is to say, the rela- 
tive population of baryons and mesons can be attributed 
not only to flow but rather to the dynamical properties 
of quark clustering in a varying density scenario. A nat- 
ural question is whether those features remain true for 
baryons and mesons with one constituent heavy-quark 
and whether a computed, as opposed to assumed, baryon 
to meson ratio, can at least partially explain the anoma- 
lous single, non-photonic electron suppression at RHIC. 

The invariant transverse momentum distribution of a 
given hadron can be written as an integral over the freeze- 
out space-time hypersurfacc E, of the relativistically in- 
variant phase space particle density F{x, P), 
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where P is the hadron's momentum, u{x) is a future 
oriented unit four-vector normal to S and g is the de- 
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generacy factor for the hadron which takes care of the 
spin degree of freedom. The function F(x, P) contains 
the information on the probabihty that the given hadron 
is formed. 

To allow for a dynamical recombination scenario in a 
thermal environment, let us assume that the phase space 
particle density F{x, P) can be factorizcd into the prod- 
uct of a term containing the thermal occupation number, 
including the effects of a possible flow velocity, and an- 
other term containing the system energy density e driven 
probability P(e), for the coalescence of partons into a 
given hadron. We thus write 

F{x,P)^e~P-'"^^^/^P{e), (8) 

where v{x) is the flow velocity. As we will show, the 
probability 7'(e) incorporates in a simple manner the in- 
formation that the coalescing partons need to be close in 
configuration space as well as to have a not so different 
velocity. 

To compute the probability 7^(e), it has been shown in 
Ref. [l^ (where we refer the reader to for details) that 
use can be made of the string flip model 0, [2l|, [2^ in or- 
der to get information about the likelihood of clustering 
of constituent quarks to form hadrons from an effective 
quark-quark interaction. In short, the model is a varia- 
tional quantum Monte Carlo simulation that, taking a set 
of equal number of all color quarks and antiquarks at a 
given density, computes the optimal configuration of col- 
orless clusters (baryons or mesons) by minimizing the po- 
tential energy of the system. At low densities, the model 
describes the system of quarks as isolated hadrons while 
at high densities, this system becomes a free Fermi gas of 
quarks. The interaction between quarks is pair-wise and 
taken as harmonic. The optimal clustering is achieved 
by finding the optimal pairing between two given sets of 
quarks of different color for all possible color charges. 

Consider, for example, two sets with equal amount A of 
quarks, one of color ci and the other of color C2, irrespec- 
tive of flavor, in accordance to the the flavor-blindness 
nature of QCD. We define the optimal pairing between 
ci and C2 quarks as the one producing the minimum in 
the potential energy, over all possible permutations: 

A 

Vc,c2 =min^v[nci,P{nc2)], (9) 

1=1 

where is the spatial coordinate of the i-th Ci quark 
and P^Tic^) is the coordinate of the mapped j-th C2 
quark. The harmonic interaction between pairs is written 
as 

•"{nc, , rjc2 ) = ^k{r,c, - Tjc^f , (10) 

where k is the spring constant. There are two possible 
kinds of hadrons that can be formed: 

i) Meson-like. In this case the pairing is imposed to be 
between color and anticolors, and the many-body poten- 
tial of the system made up of mesons is given by: 

Vmes = Vbb + VgG + VrR (H) 



where the individual terms are given by Eq. ([9|) for the 
corresponding colors. R{R), B{B) and G{G) are the la- 
bels for red, blue and green color (anticolor) respectively. 
Note that this potential can only build pairs. 

ii) Baryon-like. In this case the pairing is imposed to 
be between the different colors in all the possible combi- 
nations. In this manner, the many-body potential is: 

Vbar - Vrb + VbG + Vrg (12) 

which can build colorless clusters by linking 3(RBG), 
6(RBGRBG),... etc., quarks. Since the interaction is 
pair-wise, the 3-quark clusters are of the delta (triangu- 
lar) shape. 

According to QCD phenomenology, the formed 
hadrons should interact weakly due to the short-range 
nature of the hadron-hadron interaction. This is partially 
accomplished by the possibility of a quark flipping from 
one cluster to another. At high energy density, asymp- 
totic freedom demands that quarks must interact weakly. 
This behavior is obtained once the average inter-quark 
separation is smaller than the typical confining scale. 

To describe the evolution of a system of N quarks as a 
function of the particle density we consider the quarks 
moving in a three-dimensional box, whose sides have 
length L, and the system described by a variational wave 
function of the form: 

*a(xi,...,xw) = e-^^(''i--''")$i.G(xi,...,xjv), (13) 

where A is the single variational parameter, V^(xi,...,XAr) 
is the many-body potential defined in Eqs. pT|) and 
for mesons and baryons respectively, and $i7G'(xi,...,XAr) 
is the Fermi-gas wave function given by a product of 
Slater determinants, one for each color-flavor combina- 
tion of quarks, which are built up of single-particle wave 
functions describing a free particle in a box p^ . The 
square of the variational wave function is the weighting 
probability in the sampling, which we carry out using 
metropolis algorithm. 

The variational wave function is taken to have the form 
given in Eq. p3|) since we arc interested in the evolution 
of the system from low to high energy densities. The 
exponential term is responsible of the clustering correla- 
tions. At low energy density, the system is formed by 
isolated color-singlet hadrons and quarks strongly inter- 
acting inside each cluster; in this case, the exponential 
term of the wave function has a big contribution since 
the average interquark distance is of the order of the con- 
fining scale. In contrast, at high energy density, where 
asymptotic freedom takes place, the interaction between 
quarks is weak and the system looks like a Fermi gas of 
quarks. In this case, the inter-quark separation is much 
smaller than the confining scale and the effect of the ex- 
ponential term vanishes. Notice that these features allow 
us to identify the value of the variational parameter A as 
being directly proportional to the probability to form a 
cluster. This fact will be latter exploited to define the 
density dependent probability V{e) since, as we show be- 
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low, A changes from a fixed value at low density (isolated 
clusters) to zero at high density (Fermi gas). 



III. PROBABILITIES 

All the results we present here come from simulations 
done with 384 particles, 192 quarks and 192 antiquaries, 
corresponding to having 32 light quarks and 32 heavy 
quarks, plus their antiquarks in the three color charges 
(anti-charges). Hereafter we refer to light quarks as u- 
quarks and to the heavy ones as c-quarks. To take into 
account the mass difference between the u and c quarks 
we set rric = 10to„. We have checked that variations of 
this particular choice do not affect our relative probabil- 
ities. 

To determine the variational parameter as a function 
of density we first select the value of the particle density p 
in the box, which, for a fixed number of particles, means 
adjusting the box size. Then, we compute the energy 
of the system as a function of the variational parameter 
using the Monte Carlo method described in the previous 
section. The minimum of the energy determines the opti- 
mal variational parameter. We repeat the procedure for 
a set of values of the particle densities in the region of 
interest. To get a measure of the probability to form a 
cluster, we take the variational parameter and divide it 
by its corresponding value at the lowest density. Notice 
that since the heavy quarks are not as abundant as the 
light ones, they do not contribute to the energy density 
and thus, within the model, this last can be computed 
by assuming that only light flavors contribute. 

The information contained in the variational parame- 
ter is global in the sense that it only gives an approximate 
idea about the average size of the inter-particle distance 
at a given density, which is not necessarily the same for 
quarks in a single cluster. To correct for this, and in 
order to find an appropriate measure of the probability 
to form baryons and mesons, we need to multiply these 
variational parameters by the likelihood to find clusters 
of baryons made up of two-light, one- heavy quark and 
mesons made up of one-light, one-heavy quark. This like- 
lihood has to consider the fact that the thermal plasma is 
mainly made up of light quarks and thus that the number 
of produced heavy quarks is relatively small. To accom- 
plish this, notice that in a model where the interaction 
between quarks to form clusters is flavor (as well as color) 
blind, this likelihood should account only for the combi- 
natorial probabilities. 

Consider the case where one starts out with a set of 
n u-quarks and m c-quarks each coming in three colors. 
The number of possible colorless baryons containing three 
quarks of all possible flavors that can be formed are 

kind number 
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(14) 
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ccc 



and the total number of possible baryons is (n+m)^. The 
same counting applies for antibaryons when one starts 
from the same numbers of antiquarks instead of quarks. 

Now, consider the case where one starts with a set 
of n w-quarks, n w-antiquarks, m c-quarks and m c- 
antiquarks, each coming in three colors. The number 
of possible colorless mesons containing quark-antiquark 
pairs of all possible flavors that can be formed are 
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(15) 



and the total number of possible mesons is 3(n 4- m)^. 
We now ask for the relative abundance of baryons with 
respect to mesons computed under the above assump- 
tions on the number of light and heavy quarks that we 
start from. Since in the case of mesons we are allowing 
to consider the case uc as well as uc, wc need to include 
in the counting of the groups of three quarks also the 
antibaryons. Thus the relative abundance is 



baryons - 



antibaryons 2 x 3n^m/(ri + m)^ 
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2 X nm/ {n + to)^ 
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Let us now impose that the number of u-quarks be a 
multiple I of the number of c-quarks, namely, n = Im. 
Therefore the above relative abundance can be written 



c — baryons + c — antibaryons il 
c ~ mesons -\- c — antimesons 2(/ + 1) 



(17) 



Notice that in the plasma, the number of u-quarks greatly 
exceeds the number of c-quarks. Therefore a good ana- 
lytical estimate of the above relative abundance can be 
obtained by taking I oo which gives 
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It can be checked that the asymptotic value 3/2 is rapidly 
reached, for instance, by taking I — 30, the above fraction 
already becomes 1.475. 

Figure [5] shows the probability parameter V^'^{e) 
for baryons and mesons, obtained by multiplying the 
variational parameter with the corresponding fraction of 
baryons/mesons formed at the given energy density. In 
the case of mesons it corresponds to 1/4 irrespective of 
the density, while for baryons it has a functional form, 
since the kind of clusters can be different as density in- 
creases. For low densities the ratio of the probabilities be- 
comes 3/2, as expected from the combinatorial described 
above. Shown in the figure is also a fit to the variational 
parameters with the functional form 



f{x) = ai + 
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1 + exp [(x — x^j/dx] 
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FIG. 2; (Color online) Probabilities V^''^' to produce 
charmed baryons and mesons as a function to the energy den- 
sity e. Shown are the results of the Monte Carlo simulation 
for baryons (full circles) and mesons (open circles) together 
with a fit to these. 
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FIG. 3: (Color online) Charmed baryon and meson transverse 
momentum distributions. The parameters used in the calcu- 
lation are = 2.29 GeV, m"' = 1.87 GeV, ro = 1 fm, 
To = 200 MeV, T/ = 100 MeV, corresponding to a final time 
Tf = 8 fm. Shown are the cases with vt = and vt = 0.4. 
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We will use this analytical expression to carry out the 
calculation of the spectra that we proceed to describe. 



IV. BARYON TO MESON RATIO 

In order to quantify how the different probabilities 
to produce sets of three quarks as compared to sets of 
two quarks affect the particle's yields as the energy den- 
sity changes during hadronization, we need to resort to 
a model for the space-time evolution of the collision. 
We take Bjorken's scenario which incorporates the fact 
that initially, expansion is longitudinal, that is, along the 
beam direction which we take as the z axis and include 



transverse flow as a small effect on top of the longitudi- 
nal expansion. In this scenario, the relation between the 
temperature T and the l-fl proper-time r is given by 



T = To ('^ 



(22) 



where r = y/t'^ — z^. Equation assumes that the 
speed of sound Vs changes slowly with temperature. For 
simplicity we take Vs as a constant equal to the ideal gas 
limit ^2 ^ 1/3. 

We also consider that hadronization takes place on hy- 
persurfaces S characterized by a constant value of r and 
therefore 



dTi — Tp dp d(j) drj, 



where 
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is the spatial rapidity and p, (j) are the polar transverse 
coordinates. Thus, the transverse spectrum for a hadron 
species H is given as the average over the hadronization 
interval of the right hand-side of Eq. ([7]) , namely 
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FIG. 4: (Color online) Charmed baryon to meson ratio, Ca, 
as a function of transverse momentum. The parameters used 
in the calculation are = 2.29 GeV, = 1.87 GeV, 

7-0 = 1 fm. To = 200 MeV, T/ = 100 MeV, corresponding to 
a final time r/ = 8 fm. Shown is a range when varying the 
transverse expansion velocity vt from (upper curve at low 
momenta) to 0.4 (lower curve at low momenta). 



To find the relation between the energy density e -that 
the probabiUty V depends upon- and T, we resort to 
lattice simulations. For the case of two flavors (since the 
heavy quark does not thermalize) , a fair representation 
of the data [l7j is given by the analytic expression 



1 + tanh 



(26) 



with a = 4.82 and b = 0.132. We take = 175 MeV. 
The flow four-velocity vector is given by 



(cosh r] cosh rjT , sinh rix cos <; 
sinh riT sin 0, sinh rj cosh t^t) 



(27) 



where the magnitude of the transverse flow velocity vt 
and rjT are related by vt = tanhrjT- The normal to the 
freeze-out hypersurfaces of constant t, u'', is given by 

= (cosh ry, 0,0, sinh r/). (28) 

We write the momentum four-vector in components as 

P'^ = (my cosh y,pT cos $,pt sin $, TOj- sinh y), (29) 

where y is the 1 + 1 momentum rapidity given by 
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and <I> the azimuthal angle of the momentum components 
in the transverse plane. 

Therefore, the products P ■ u and P ■ v appearing in 
Eq. can be written as 

P ■ V = rriT cosh(77 — y) cos rjT — pt cos(^ — <i>) sinh tjt 
P-u = niT cosh{r] — y) , (31) 

Considering the situation of central collisions, we can as- 
sume that there is no dependence of the particle yield on 
the transverse polar coordinates. Integration over these 
variables gives 
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(32) 



where p„uci is the radius of the colliding nuclei and Iq is 
the Bessel function / of order zero. 

We now consider as a further simplification that the 
space-time and momentum rapidities are completely cor- 
related, that is rj ^ y. Under this assumption, the inte- 
gral over T] in Eq. (j32|) can be performed and we finally 
get 
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Armed with the expression to compute the hadron trans- 
verse momentum distribution, wc now proceed to apply 
the analysis to the computation of the charmed meson 
and baryon distributions. 



V. RESULTS 

Figure [3] shows examples of the transverse momen- 
tum distributions for mesons and baryons obtained from 
Eq. ([55)1 . We set the masses of the charmed baryons 
and mesons as = 2.29 GeV (corresponding to Ac) 
and m^^ = 1.87 GeV (corresponding to D). We take 
the initial hadronization time as tq = 1 fm, at an initial 
temperature To = 200 MeV and the final hadronization 
temperature as Tj = 100 MeV, corresponding, according 
to Eq. ((22)) . to a final time tj = 8 fm. Shown are the 
cases with vt = and vt — 0.4. Notice that a finite 
transverse expansion velocity produces a broadening of 
the distributions, as expected. 

Figure 3] shows the charmed baryon to meson ratio ob- 
tained from the ratio of the above transverse momentum 
distributions. Shown is a range for this ratio when vary- 
ing the transverse expansion velocity vt from to 0.4. 
Notice that for a finite vt, this ratio goes above 1 for 
PT > 3.5 GeV. 
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FIG. 5: (Color online) Suppression factor T%j^ as a function 
of transverse momentum. The parameters used in the cal- 
culation are = 2.29 GeV, ml^ = 1.87 GeV, ro = 1 fm, 
To = 200 MeV, T/ = 100 MeV, corresponding to a final time 
T/ = 8 fm, X = 0.14, a=0.073. Shown is a range for the trans- 
verse expansion velocity form = (upper curve at low pr) 
and vt ~ 0.4 (lower curve at low pr- 



We now proceed to compute the pr unintegrated func- 
tion T^^. For this purpose, we take that the pos- 
sible charmed mesons decaying inclusively into elec- 
trons or positrons arc D± ^ 16.0%), D°,Z)° 



6.53%), Df {B^ 



and that 

the possible charmed baryons decaying inclusively into 



electrons or positrons are Ac, Ac {B 
Thus we get 

X = 0.14. 



4.5%). 



(34) 



We also approximate the masses of all the charmed 
mesons considered to be equal to the mass of the 
mesons. 

From Eq. ^ we see that, without integrating over 
Pt, the dependence on the transverse momentum comes 
from a = {dN^^J dpx) I {dN^^J dpx) and the product Ca = 
{dN'UdpT)l{dN^JdpT\_^^c integrated ratio a'"' has 
been computed in Ref. [l^l using a Pythia simulation, 
with the result a'"' = 0.073. We have also performed 
a simulation using Pythia at NLO with 100,000 events 
and have found that with such statistics, the ratio of 
charmed baryons to charmed mesons in p -I- p collisions 
at y/sNN = 200 GeV is flat up to — 5 GeV and con- 
sistent with the value reported in Ref. [l3- Therefore, 
for simplicity we take a as a constant equal to the above 



(1 + xa'-) 1 + {dN^JdpT)/{dN^JdpT) 



Figure O shows as a function of pr ■ We have used 
a range of values for the transverse expansion velocity 
between wt = and vt = 0.4. We see that for the 
chosen evolution parameters, is indeed smaller than 
1 and thus it contributes to the suppression of the single 
non- photonic electron nuclear modification factor i?^^. 



VI. CONCLUSIONS 

In this work we have shown that the anomalous sup- 
pression of the single non-photonic electron nuclear mod- 
ification factor R\j^ can be partially understood by re- 
alizing that this quantity is affected by an enhancement 
in the charmed baryon to meson ratio at intermediate 
Pt in Au -I- Au collisions. This enhancement is due to 
the fact that in this region, thermal recombination be- 
comes the dominant mechanism for hadron production. 
We have made use of the DQRM to calculate this ra- 
tio and have shown that for moderate and even for van- 
ishing transverse expansion velocities, it indeed can be 
larger than the charmed baryon to meson ratio in p -I- p 
collisions. This enhancement in turn produces that the 
function is below 1 and thus contributes to the sup- 
pression factor introduced by considering energy losses 
due to the propagation of heavy flavors in the plasma. 

It is worth to keep in mind some important features 
concerning the results of the present calculation: First, 
notice that we have not included the momentum shift in- 
troduced by energy losses when computing the transverse 
distributions of charmed mesons and baryons. This is so 
because for i?^^ , energy losses should be included in the 
prcfactor of the function T|^. In this sense, in order to 
avoid a double counting of the effect, the ratio that goes 
into the calculation of this last function is the raw ratio. 
Second, it is expected that at some value of pt, frag- 
mentation becomes the dominant mechanism for hadron 
production and therefore that the charmed baryon to me- 
son ratio decreases above that pT value, given that frag- 
mentation produces more mesons than baryons. Third, 
we have considered flnite values of transverse flow for 
charmed mesons and baryons even thought it might be 
questionable that heavy flavors also flow as light flavors 
do. Nevertheless, there seems to be some experimental 
support for heavy quark flow p3j . In this sense, the flow 
strength range we have considered is only for moderate 
values. Notice however that even in the absence of flow 
the suppression factor keeps being less than 1. Some of 
these issues will be the subject of a future work to appear 
elsewhere. 
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